Double-exchange theory of ferroelectric polarization in orthorhombic manganites with twofold periodic magnetic texture

We argue that many aspects of improper ferroelectric activity in orthorhombic manganites can be rationalized by considering the limit of infinite intra-atomic splitting between the majority- and minority-spin states (or the double exchange limit), which reduces the problem to the analysis of a spinless double exchange (DE) Hamiltonian. We apply this strategy to the low-energy model, derived from the first-principles calculations, and combine it with the Berry-phase theory of electric polarization. We start with the simplest two-orbital model, describing the behavior of the eg bands, and apply it to the E-type antiferromagnetic (AFM) phase, which in the DE limit effectively breaks up into one-dimensional zigzag chains. We derive an analytical expression for the electronic polarization (Pel) and explain how it depends on the orbital ordering and the energy splitting Delta between eg states. Then, we evaluate parameters of this model, starting from a more general five-orbital model for all Mn 3d bands and constructing a new downfolded model for the eg bands. From the analysis of these parameters, we conclude that the behavior of Pel in realistic manganites corresponds to the limit of large Delta. We further utilize this property in order to derive an analytical expression for Pel in a general two-fold periodic magnetic texture, based on the five-orbital model and the perturbation-theory expansion for the Wannier functions in the first order of 1/Delta. This expression explains the functional dependence of Pel on the relative directions of spins. Furthermore, it suggests that Pel is related to the asymmetry of the transfer integrals, which should simultaneously have symmetric and antisymmetric components. Finally, we explain how the polarization can be switched between orthorhombic directions a and c by inverting the zigzag AFM texture in every second ab plane.Comment: 41 page, 10 figure

Spin dependence of ferroelectric polarization in the double exchange model for manganites

The double exchange (DE) model is systematically applied for studying the coupling between ferroelectric (FE) and magnetic orders in several prototypical types of multiferroic manganites. The model was constructed for the magnetically active Mn $3d$ bands in the basis of Wannier functions and include the effect of screened on-site Coulomb interactions. The essence of our approach for the FE polarization is to use the Berry phase theory, formulated in terms of occupied Wannier functions, and to evaluate the asymmetric spin-dependent change of these functions in the framework of the DE model. This enables us to quantify the effect of the magnetic symmetry breaking and derive several useful expressions for the electronic polarization ${\bf P}$, depending on the relative directions of spins. The proposed theory is applied to the solution of three major problems: (i) The magnetic-state dependence of ${\bf P}$ in hexagonal manganites; (ii) The microscopic relationship between canted ferromagnetism and ${\bf P}$ in monoclinic BiMnO$_3$; (iii) The origin of FE activity in orthorhombic manganites. We show that for an arbitrary noncollinear magnetic structure, propagating along the orthorhombic $\boldsymbol{b}$ axis and antiferromagnetically coupled $\boldsymbol{c}$, ${\bf P}$ can be obtained by scaling the one of the E-phase with the prefactor depending only on the relative directions of spins and being the measure of the spin inhomogeneity. This picture works equally well for the twofold (HoMnO$_3$) and fourfold (TbMnO$_3$) periodic manganites. The basic difference is that the twofold periodic magnetic structure is strongly inhomogeneous, that leads to large ${\bf P}$. On the contrary, the fourfold periodic magnetic structure can be viewed as a moderately distorted homogeneous spin spiral, which corresponds to weaker ${\bf P}$.Comment: 32 pages, 7 figure

Microscopic analysis of the magnetic form factor in low-dimensional cuprates

We analyze the magnetic form factor of Cu$^{2+}$ in low-dimensional quantum magnets by taking the metal-ligand hybridization into account explicitly. In this analysis we use the form of magnetic Wannier orbitals, derived from the first-principles calculations, and identify the contributions of different atomic sites. Having performed local density approximation calculations for cuprates with different types of ligand atoms, we discuss the influence of the on-site Coulomb correlations on the structure of the magnetic orbital. The typical composition of Wannier functions for copper oxides, chlorides and bromides is defined and related to features of the magnetic form factor. We propose easy-to-use approximations of the partial orbital contributions to the magnetic form factor in order to give a microscopic explanation for the results obtained in previous first-principles studies.Comment: 5 pages, 4 figure

Validity and limitations of the superexchange model for the magnetic properties of Sr2IrO4 and Ba2IrO4 mediated by the strong spin-orbit coupling

Layered perovskites Sr2IrO4 and Ba2IrO4 are regarded as the key materials for understanding the properties of magnetic relativistic insulators, mediated by the strong spin-orbit (SO) coupling. One of the most fundamental issues is to which extent these properties can be described by the superexchange (SE) model, formulated in the limit of the large Coulomb repulsion. In the present work we address this issue by deriving the relevant models and extracting parameters of these models from the first-principles calculations. First, we construct the effective Hubbard-type model for the t2g bands, by recasting the problem in the language of Wannier orbitals. Then, we map the obtained electron model onto the pseudospin model by applying the theory of SE interactions. We discuss the microscopic origin of anisotropic SE interactions, inherent to the compass Heisenberg model, and the appearance of the antisymmetric Dzyaloshinskii-Moriya term, associated with the additional rotation of the IrO6 octahedra in Sr2IrO4. In order to evaluate the Neel temperature (TN), we employ the non-linear sigma model. While for Sr2IrO4 our value of TN agrees with the experimental one, for Ba2IrO4 it is overestimated by a factor two. We argue that this discrepancy is related to limitations of the SE model: while for more localized t2g states in Sr2IrO4 it works reasonably well, the higher-order terms, beyond the SE model, play a more important role in the more "itinerant" Ba2IrO4, giving rise to the new type of isotropic and anisotropic exchange interactions. This conclusion is supported by unrestricted Hartree-Fock calculations for the same electron model, where in the case of Ba2IrO4, already on the mean-field level, we were able to reproduce the experimentally observed magnetic ground state, while for Sr2IrO4 the main results are essentially the same as in the SE model.Comment: 37 pages, 9 figure

Lattice Distortion and Magnetism of 3d-t2gt_{2g} Perovskite Oxides

Several puzzling aspects of interplay of the experimental lattice distortion and the the magnetic properties of four narrow $t_{2g}$-band perovskite oxides (YTiO$_3$, LaTiO$_3$, YVO$_3$, and LaVO$_3$) are clarified using results of first-principles electronic structure calculations. First, we derive parameters of the effective Hubbard-type Hamiltonian for the isolated $t_{2g}$ bands using newly developed downfolding method for the kinetic-energy part and a hybrid approach, based on the combination of the random-phase approximation and the constraint local-density approximation, for the screened Coulomb interaction part. Then, we solve the obtained Hamiltonian using a number of techniques, including the mean-field Hartree-Fock (HF) approximation, the second-order perturbation theory for the correlation energy, and a variational superexchange theory. Even though the crystal-field splitting is not particularly large to quench the orbital degrees of freedom, the crystal distortion imposes a severe constraint on the form of the possible orbital states, which favor the formation of the experimentally observed magnetic structures in YTiO$_3$, YVO_, and LaVO$_3$ even at the HF level. Beyond the HF approximation, the correlations effects systematically improve the agreement with the experimental data. Using the same type of approximations we could not reproduce the correct magnetic ground state of LaTiO$_3$. However, we expect that the situation may change by systematically improving the level of approximations for dealing with the correlation effects.Comment: 30 pages, 17 figures, 8 tables, high-quality figures are available via e-mai

Stability of an oscillating boundary layer

Levchenko and Solov'ev (1972, 1974) have developed a stability theory for space periodic flows, assuming that the Floquet theory is applicable to partial differential equations. In the present paper, this approach is extended to unsteady periodic flows. A complete unsteady formulation of the stability problem is obtained, and the stability characteristics over an oscillating period are determined from the solution of the problem. Calculations carried out for an oscillating incompressible boundary layer on a plate showed that the boundary layer flow may be regarded as a locally parallel flow